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Abstract 

Motivated by recent results on the index of the Dirac operator D = "j^D^ of QCD on lattice 



and also by results on topological features of electrons and holes of 2-dimensional graphene, we 
compute in this paper the Index of D for fermions living on a family of even dimensional lattices 
denoted as 'L2N and describing the 2N-dimensional generalization of the graphene honeycomb. The 
calculation of this topological Index is done by using the direct method based on solving explicitly 
the gauged Dirac equation and also by using specific properties of the lattices 'L2N which are shown 
to be intimately linked with the weight lattices of SU {2N + 1). The Index associated with the 
two leading = 1 and N = 2 elements of this family describe precisely the chiral anomalies of 
graphene and QCD4. Comments on the method using the spectral flow approach as well as the 
computation of the topological charges on 2-cycles of 2N-dimensional compact supercell in L,2N 
and applications to QCD/j^ are also given. 

Key words: Topological index, Lattice QCD4, Graphene, Chiral anomaly. Root and Weight 
lattices of SU (N), Adam's spectral flow. 
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I. INTRODUCTION 



It is quite well known that the spectrum of the gauged Dirac operator D = 7^ {d^ — iA^) 
= ^^D^ of 2- and 4- dimensional Dirac theories describing the dynamics of fermionic waves 
\1/ (x) in a uniform background field with topological charge Qtot obeys the famous Atiyah- 
Singer index theorem [l], I2I; one of the most substantial achievements of modern mathemat- 
ics. The index of the Dirac operator D, to which we refer below to as Ind{D), relates the 
topological charge Qtot to the net numbers and of chiral and antichiral zero modes 
of D as 

Qtot = N+-N_ (1) 

showing in turns that the background field breaks implicitly the left-right parity symme- 
xy. The Ind {D) is a powerful topological quantity that has been used for diverse purposes 



3; 



see also 



15[ and refs therein for other applications; it gives a rigorous explanation 



m 

of the origin of chiral anomalies and constitutes an alternative approach to the perturba- 
tive method based on computing radiative corrections to describe ^'-f^'^Afj^ interactions in 
QFTi^2 where the topological Chern-Simons gauge theory emerges as a 1-loop correction 
of the gauge field propagator [la, ll7|. In 2-dimensions, the computation of Ind{D) shows 



that the anomalous quantum Hall effect (QHE) of graphene |l8l42Cll| is precisely due to the 
chiral anomaly of zero modes; a basic result that is expected to be valid as well for higher 
dimensional Dirac fermions in background fields including light quarks on 4D hyperdiamond 
2l| - 28| and fermions on higher 2N- dimensional honeycombs. Recall by the way that in 
graphene the quantum Hall effect is a very special effect in the sense it can be observed at 
room temperature] the gap between the n = and n = 1 Landau levels is around 1300K 
ai 10 Tesla compared to around 1 OOK in an ordinary 2D electronic gaz 3ll, l32| ; see also 



33 



36| for other related aspects. 



Motivated by recent developments on the index theorem on lattice QCD4 37H40|], we com- 
pute in this paper the Ind{Di), and in general the Ind{D2N), of fermions in background 
fields living on a class of even- dimensional lattices L2Ar describing the 2N- dimensional 
generalization of the honeycomb L2. The calculation of Ind{D2N) will be done as follows: 

• use known results on the index to bring the lattice analysis to the spectrum of the 
fermionic operator near the Fermi level; i.e the spectrum near the Dirac points of 
lattice QCD2Ar where live fermionic zero modes contributing to Ind{D2N)- 



• then determine these fermionic zero modes by using the direct method based on the 
two following: 

(a) using specific features of the SU {2N + 1) roots and weights to deal with the 
symmetry properties of the lattices with = 1, 2, 

{b) solving explicitly the gauged Dirac equation D"^ = e\l', with e standing for massive 
deformations around the Fermi level. 



Recall that there are two basic ways to compute Ind^D^); {{) the direct method which we 
will be considering here; and (ii) the so called spectral flow approach of Adams 37H40| 
based on the introduction of a hermitian version Hgp of the Dirac operator. In QCD4 on 
hyperdiamond, this spectral hamiltonian has the form Hgp = 75 (D — m) showing that any 
zero mode of D with ± chirality corresponds to some eigen-modes of Hsp with eigenvalues 



A (m) = ^m; for details see [28|, l37H42|. 



Our interest into the study of the index of the Dirac operator of fermions on has been 
also motivated by the two following: 



(1) the 2 leading lattices L2, L4 of the family are precisely given by the honeycomb of 
graphene, thought of here as lattice QCD2, and the 4D hyperdiamond used in lattice QCD4. 
So, one expects the members of this family to share some basic features; in particular 
methods to approach the physical properties of fermions on L,2n- This link opens a window 
for new ways to modeling and simulatin g in euclidian relativistic theory by borrowing ideas 



from graphene as done by M. Creutz in 2l|]; see also 22 



23|, 



25|. 



(2) the existence of a remarkable relation between the IL2N honeycomb's class and the fam- 
ily of weight lattices of the SU {2N + 1) Lie algebras. This extraordinary relation allows 
a unified description of the tight binding description of fermions on h2N with a generic 
integer; and permits moreover to take advantage of the power of the SU {2N + 1) repre- 
sentations to work out explicit configurations for fermions and gauge fields on h2N- For 
example boundary conditions of the fields on supercells in h2N as well as the determination 
of the Dirac points {Pi} regarding fermions on h2N get mapped to manageable equations on 
the SU {2N + 1) weight lattice where the duality between siniple roots CKj of SU {2N + 1) 



and its fundamental weight vectors Ui plays a crucial role 



271,1291,13 



Our explicit analysis for computing the index of the Dirac operator allows us as well to get 
more insight into the structure of the topological index on lattice; in particular into the two 
following things: (a) the relation between the various possible fluxes Qi through the 2-cycles 
Ci of the supercell compactification and the total charge Qtot of the background fields, (b) 
the role played by the different matrices '^2n+i construct from the 2D Pauli ones 

namely 

T,,^,=T,(^{a,r^®...(^{a,r^-^ (2) 
with nj = 0, 1 and Ts = diag (1, —1) = as. 

Regarding the topological charges Qi and Qtot, notice that the lattice 'L2N is recovered by 
considering a 2N- dimensional compact supercell with homology classes as those of the real 
2N-torus T^^. While in 2-dimensions the homology of the supercell has one 2-cycle given 
by the parallelogram of fig El the situation is richer in higher dimensions. In the case of 4D 
hyperdiamond for instance, the homology of the 4D compact supercell has, in addition to 
the real 4-cycle C4 ~ T"^, a basis of six 2-cycle Cj ~ leading to the following topological 
charges 

related as in eg dlSSD and playing a central role in computing the index of the Dirac operator. 

Concerning the matrices ^^n+i fermions on 2N-honeycombs, one has in general 2^~^ 
possible representations depending on the values of the n^'s of ([2]). For the case of the 4D 
hyperdiamond, there are two kinds of such matrices and are as follows 

75 = 0-3®/, 0-3 ®r3. 

The explicit computation of the Dirac index given in sections 5 and 6 shows that these 
matrices lead to different relations between the zero modes and the topological charges. 
We will show that the right index of the Dirac operator that recovers ([1]) is given by 
Tr {^a^ ® Ta^) = Q^p- 

The presentation is as follows: In section 2, we review the usual approach to honeycomb; 
but now by using the root and weight lattices of SU (3); the latter is a hidden symmetry of 
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the honeycomb. In section 3, we study the topological aspects of fermions in graphene in 
presence of an external magnetic field B and develop the explicit computation of the zero 
modes of the gauged Dirac operator. As we will show, the basic properties of the ground 
state are encoded in the sign of the background field. In section 4, we use the results on 
graphene to approach the fermions on hyperdiamond in presence of two background fields 
B and E. In this study, we also use roots and weight lattices of SU (5) that appears as a 
hidden symmetry of the 4D crystal. In section 5, we compute the spectrum of the Dirac 
operator in continuum, and in section 6 we work out the complete spectrum of fermions 
on the hyperdiamond lattice. We also compute the topological index giving the relation 
between the zero modes and the various topological charges of the background fields. Here 
also, we show that the ground state features are encoded in the sign of and -j^. In section 
7, we give a conclusion and make comments on the extension to fermions on 2N- dimensional 
honeycomb in presence of N background fields Bi whose signs characterize completely 
the ground states. 



Before proceeding, we would like to notice that we will omit certain technical details on 
lattice calculations from the core of the paper; the essential of these details is reported in 



appendix 8; and related extensions can be found in 21 



-25 
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II. FERMIONS ON 2D HONEYCOMB 



In this section, we develop the so called primitive compactification of fermion on honey- 
comb used in 20| to study topological aspects of fermions in graphene. To get more insight 



into our analysis, we first describe the canonical frame formulation, generally used to deal 
with fermions on square lattice; then we turn back to study the primitive frame (P-frame) 
and exhibit its link with the root and the weight lattices of the SU (3) Lie algebra. This 
study can be also viewed as a step to fix the ideas before moving to the case of QCD4 
fermions on hyperdiamond to be considered in sections 4, 5 and 6. 



A. 2D honeycomb and SU (3) weight lattice 

Honeycomb is a 2D lattice where each site has 3 first nearest neighbors and 6 second 
nearest ones as shown on fig [H the 3 first neighbors will be viewed here as the vertices of 
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an equilateral triangle transforming as a SU (3) triplet; and the second 6 ones the vertices 
of an hexagon. The lattice parametrization of the honeycomb can be obtained from the real 
plane by restricting the local real coordinate r = {X,Y) to r„,„ = (X„, Fm) with n, m 
integers. As a 2D vector, one ma y ex pand the position vector r in diverse ways; in particular 
into 2 particular frames used in 20| and named there as perpendicular frame and primitive 
one; see fig [1] for illustration. Here, we will refer to the perpendicular frame as the canonical 



one. 





FIG. 1: On left the perpendicular compactification. On right the primitive compactification 



• Canonical compactification 

In the canonical compactification (C-frame), one uses the usual orthogonal cartesian basis 
{ex, ey} of to decompose vectors like r = Xex + Yey, 

e^.e^ = , X = ex.r , Y = ey.r . (3) 

The restriction of local fields F (r) to the honeycomb is obtained by using the Dirac delta 
function S {r — rn,m) as 

i^r„.„. = / dh5ir- r„,^) F (r) . (4) 

The dynamics of fermionic waves \1/ (T, X, Y) describing the delocalized electrons near the 
Dirac points is given by the Dirac equation 



with 

satisfying the ChfFord algebra in (2 + l)-dimensions. The solutions of this equation are 
given by plane waves e^'^^'^^'^'^^^+'^y^) with = + ky- Boundary conditions such as 
^ {T,X + L,Y) = (T, X, Y) leads to discrete Kx = ^,neZ. 

• Primitive compactification 

In the primitive frame (P-frame) of honeycomb, one uses the non orthogonal vector basis 

allowing to take advantage of the two following remarkable features of the honeycomb: 

(i) honeycomb contains the root lattice of SU (3) 
The basis vectors e^,, e^^ are, up to a scale factor, precisely the simple roots ai and 0:2 of the 
SU (3) Lie algebra 

so that space vectors like positions can be decomposed as 



r = x^ai + x^a^ ■ (8) 
The relation between the C- and P- frames is given by 



or by using matrices 



The inverse transformation is 



V2 ( al aj 
al a: 



(9) 





(10) 





with the constraint /^^a^'' = 5^ which is solved in terms of the af's hke 



with A = det (af ) 

(ii) honeycomb is the weight lattice of SU (3) 
The matrix is precisely the entries of the fundamental weight vectors of SU (3) 



These vectors are dual to the simple roots 



a,- .w 



' -SI 



Notice that the roots of Lie algebras can be expanded in terms of the weight vectors and 
vice versa. For the example of SU(3), the simple roots decompose as 



This features teaches as that honeycomb sites are given by riiuj'^ + ?7,2W^ with rii arbitrary 
integers. Moreover sites in the sublattices A and B of honeycomb are rather given by 
riiai + n2a2. 

(iii) from 2D honeycomb to 4D hyperdiamond 
The link between honeycomb and the root lattice of SU (3) is very suggestive. It allows to 
extend results on 2D graphene to fermions on higher 2N- dimensional lattices. The latters 
are isomorphic to the weight lattice of SU {2N + 1) Lie algebras and permits the following 
correspondence 



N 


SU {2N+1) 


Lattice 


1 


SU (3) 


2D honeycomb of graphene 


2 


SU{5) 


4D honeycomb = 4D hyperdiamond of QCD4 


3 


SU{7) 


6D honeycomb 



In this view, 2D graphene appears as the leading term of a family of lattice models. The 
second element of this family is remarkably given by the 4D hyperdiamond that is used in 
dealing with quarks in lattice QCD4 [2l|-[28|. 



To make contact with the study of [20|, we take the components of the simple roots in the 
C-frame as, 

ai=(^,f),«2 = (-v^,0), (12) 
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leading to 



2 2 

-V2 



/3, 



with «! A ^2 = A = -\/3; and 



a; 



a; 





" 2 
3 6 



V2 VE 

2 ' 6 



(13) 



(14) 



with ui Au2 = ^ which is | smaller compared to ai A 02. Notice that the transpose vectors 



a; 



[at 



,b2 = (/3^) = 



f,0 



3 ' 6 



(15) 



satisfy as well the duahty duality relation a^.b^ = 6^. Notice also that the vectors 0, 1, 2 
used in [2(| are nothing but the three 3 roots SU(3) namely ai, 02, = ai + a^^ see also 
fig [2] for illustration. The area of the unit cell in the C-frame is given by -L1L2 and in the 



A 


A 




B 




unit cell 




FIG. 2: The real unit cell of the 2D honeycomb in C- and P- frames. Sites A (resp B) belong to 
sublattice A (resp B) of the honeycomb. 



P-frame is hh^- The equality of these areas can checked by using eqs fl9][T2|) 



B. Dirac equation in primitive frame 

One of the lessons we have learned from above is that simple roots and fundamental 
weight vectors of are appropriate tools to deal with fermionic waves on honeycomb and 
reciprocal space. Real positions r and wave vectors k in the momentum space can be 
decomposed as, 

r = x^ai , X* = u\r , , 

(16) 

k = kiu"^ , ki = ftj.k 
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The respective norms read as = Cijx^x^, = G^^kikj with metrics Qj = aj-cuj and 
G^^ — co^.u!^ given by 

c^J =1 \ ,G-=\ yt \ (17) 




-1 2 

Plane waves e**^"" which read in C-frame as e''^'^^'^ takes also the form e**^'^' in the P- 
frame because of the duality relation oj\aj — Sj appearing in expi [uj\aj kiX^. The Fourier 
transform reads as usual 

/(r) = / (18) 

and the gradients are related as — and -A- = -jjl-. The two sublattices A and 

B of the 2D honeycomb are as follows 

A : r„ = niSii + n2a2 

B : r„ = (m - \) ai + (na - |) 

where = d^J^ai = -^oii with (i standing for the lattice parameter and n — irix^n-i) e 1? . 
The area of the unit cells in real and momentum spaces are given by 



(19) 

3(F 3 ^ 3 



laiAaal =^ = a , 

lU aK I — V3 _ 87r2 V3 
D1AD2 — 0J2--0- — -;jr-o- 



with bj = Icj' = ^^^ci^' the dual vectors of a^. For a real supercell 

SC2 = (iViai) X (A^2a2) 



described by a parallelogram with edges Ij = A^ja^, the real area S is given by <S = A^iTVa^^-^- 
The dynamics of fermionic 
it reads in the P-frame as, 



4 

The dynamics of fermionic waves in the C-frame is given by the Dirac equation 7^^^^^ = 0; 



^°|- + ^'ir- + ^'|-)* = o (20) 

at ox ay J 
with the new Dirac matrices related to the 7'*'s hke 

= uji^Y , Y = , TO = 70 (21) 
and obeying the following Chfford algebra 

{T\ T^} = 2G'^ , {TO, T^} = 
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III. TOPOLOGICAL ASPECTS OF FERMION ON HONEYCOMB 



First we give some general features regarding fermion on supercell; then we study the 
spectrum of the Dirac equation of these fermions in a uniform background field B. After 
that, we compute the index of the Dirac operator in presence of B. 

A. Fermions on supercell 

To study the topological aspects of fermions on honeycomb, one has to transcribe the 
expressions of the usual fields \t' (x) and (x) in the continuum to the case of the lattice. 
This extension is not straightforward because the field expressions seem to be in conflict 
with the lattice periodicity. Below, we study this issue by working in the P-frame of the 
honeycomb and by using its underlying SU (3) symmetry. 

1. Supercell in honeycomb 

The sites of a honeycomb supercell SC2 are of two types: A-type and B-type; and are 
defined as, 

SC2 = {r^ = niai + n2a2 } U = (m + ^)ai + (^2 + ^)a2 | (23) 

with n — (711,712) integers restricted as 

<ni< Ni-1 , Li^ dVSNi = aNi . (24) 

The vectors and are the lattice positions along the ai and a2 directions with 

A^i, N2 the respective number of sites in these directions. So the total number of the A-type 
sites in the supercell is A^iA^2 and then the total area S of the supercell is equal to NiN2^^. 
Notice also that in the supercell SC2 we also have B-type sites r"^^ with same number as 
the A-type ones so that the total sites in the supercell is 

2N1N2. (25) 

Notice that the r®^^ sites are shifted with respect to the sites by a constant vector 

as r" = -I- s with s =|ai -|- |a2. This shift vector s has two remarkable features: first it is 
independent of the site positions (711,712); and second it is proportional to — |q;i -|- \a2] 
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FIG. 3: A supercell in primitive compactification; it has 2J^ sites: 12 sites are A-type (balls in 
red) and the 12 others of B-type (balls in blue). The origin of the primitive frame belongs to 
A-sublattice. A shift by the vector s maps A-type sites into B-type ones. 




we have s = d^'^u^ . Using the parametrization of [20|, one can check that iS^ = (0, ^) 
and by substituting back into the expression of the shift vector, we get 



(26) 



2. Gauge field and fermionic waves on supercell 

To deal with the Dirac equation on honeycomb supercell in presence of the background 
field B, one has to worry about the boundary conditions of the gauge field and the fermionic 
waves under lattice periodicity 

r„ — )• r„ + nA^'iai + mA^2a2, n, m G Z (27) 



Constructions regarding this issue has been first done by Smit, J. C. Vink in [5] by working 



in the C-frame of a square lattice; see also 



20 







281]. Below we extend this analysis to the 



P-frame; this is helpful when we consider the extension to QCD4 on 4D hyperdiamond and 
extension to higher dimensional lattices. 

• Periodicity of the gauge field 

In the canonical frame, the gauge field is given by the vector while in the primitive 
frame, it is denoted as Ai. These two fields capture the same physical property; they are 
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related as = u'^^Ai or inversely Ai = ct'^A^. 

Here, we use the P-frame coordinates = {x, y) and work with the gauge covariant deriva- 
tives Vi = ^ — iAi with curvature 

[Vi,Vj] = -iTi, =ie,jB (28) 

related to the C-frame one [D^, DJ\ = —iF^^ like 

T- ■ = a^a""- F 

As usual, the gauge field Ai is defined up to gauge transformations A^ = Ai + iVLdiVL~^ 
where Vt = e"*^'^'') is a f/ (1) gauge group element with gauge parameter A (x). We also have: 

^ ^ ^ ^ (29) 

Af =Ai + dA 

Since B is constant, a typical gauge configuration of Ai that solve ( 128|) in the continuum is 
as follows, 

A (x, y) = \By , Ay (x, y) = -\Bx . (30) 

This not unique since a different, but equivalent, gauge field representation with curvature 
B is given by the following simpler one used by Smit, J. C. Vink in jsl , 

A^{x,y) = Q , Ay{x,y) = -Bx . (31) 

However, this gauge field configuration, which is valid in the continuous case, has to be 
adapted to the lattice; neither (130|1 nor ( 13T|) is periodic on the honeycomb since 



Ai {x + Li,y) ^ Ai {x,y) , 
Ai {x,y + L2) ^ Ai {x,y) . 

Explicitly, by using the expression (!30|) for the gauge field Ai, we have 

•^x (^1 y) \y=L2 2 ^ ' \y=o ' 

Ay (x, ?/) |x=Li 2 ' ^) |2;=0 • 



(32) 



(33) 



If using the gauge configuration (I3T1) . we also have 



Ax (X, y) \y = L2 =0 , A' (X, 2/) \y=0 = , 



(34) 
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To overcome this difficulty, we use the gauge symmetry freedom to ensure field's periodicity. 
This is done by requiring 

A' (X, y) \y=Q = A (X, y) \y=L2 + , .^^x 

(35J 

A y) U=o = A y) U'=Li + Qy^ , 

where now the gauge group parameter A [x, y) is no longer an arbitrary function as in the 
continuum; this is the price to pay to implement lattice periodicity. 

In what follows, we choose the gauge configuration fl3T]) used by Smit and Vink for the case 
of gauge fields on square lattice. This choice is tricky and allows tremendous simplifications 
when computing the spectrum of the Dirac operator; especially the degeneracy of the ground 
state. Using this choice the gauge covariant derivatives reduce to, 

l^x , Qy i-^y , 

with no gauge component along the x-direction. In this case, the periodicity condition on 
the gauge field namely 

Sy {X, y) \x=0 = Qy {X, y) |a;=Li + dy\ 

leads to dy\ = BLi and then to A = BLiy + Lp{x) = (y) + Lp (x) with Lp a real function in 



the variable that we shall ignore below as it doesn't affect the analysis |47|]. Notice that 
as a function on lattice, the gauge symmetry element Q, thought of as a function of the 
y-variable only, i.e: f2 = f2 (y), should be also periodic Q{y + L2) = {y). The solving of 
this condition requires the quantization of the background field as 

BL1L2 = 2itQ , Q G Z . (36) 

So the gauge group element reads as 

n (y) = exp {-t'i^y) , (37) 

and can be thought of as describing a wave plane propagating along the y-direction with the 
momentum Ky = ^j^. This property teaches us that gauge transformations are generated 
by the shifts 

Py^Py + n— — , neZ (38) 
of the momenta Py = ^ along the y-axis. 

• Periodicity of fermionic waves 
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Fermionic waves \1/ {x, y) on honeycomb supercell have to obey consistent boundary con- 
ditions. Because of gauge freedom, the boundary condition should be as 

^ (x + nLi,y + mLa) = 7^„,m^' (x, y) , 

where TZn,m is some representation of the gauge symmetry which is given by some local 
phase e^'^"-'"^^\ However, with the gauge choice (|3T|) . these condition reduces to 

'^{x,y = L2) = ^ (x, ?/ = 0) , ^^^^ 

^ (x = Li, y) = e^^'^y'^^ {x = 0,y) . 

To get the explicit expression of the fermionic waves, one has to solve the Dirac equation on 
the lattice; this will be done in next subsection; for the moment we use general arguments 
to derive some useful information on these waves. 

First, because of the choice (13T!) where the gauge field Ai has no dependence in the y- 
variable, the periodicity condition of fermionic waves namely \E' (x, y + L2) = \E' (x, y) can 
be solved in terms of plane waves propagating along the y-direction as follows 



^ (x, y) = ^e'^^^y^i (x) (40) 

The Fourier modes (x), which depend on x, carry an integer charge / and because of the 
second condition of (139!) are expected to be not completely independent fields. 
Second, to solve the boundary condition along the x-direction, we expand the wave function 
at X + Li in a similar manner as in fHOl) 



^(x + Li,y) = ^e*^^^;(x + Li) (41) 

i& 

and then require that the value of the fermionic waves at the equivalent positions x and 
X + Li to be related by a gauge transformation like \E' (x, y) = q-^^^^v v]/ (x + Li, y). This 
leads to the following relation between the Fourier modes, 

(x + U) = ^z_Q (x) (42) 

This property can be exhibited by help of the quantization relation 2'kQ = BL1L2 allowing 

i — V 

to express the wave plane basis e in terms of the magnetic field and the topological 
charge like e^^^^Q^ . Substituting this expression back into fHOj) . one ends with f H2|) . Notice 
that eq(H2l) is solved by taking 

(x) = $ (^x + ^Li^ (43) 
15 



showing that a translation y y + L2 along the y-direction by one period induces the 
shift X — )■ X + along the x-axis. The explicit expression of the function $ (^) will be 
determined below. 



B. Index of Dirac operator and chiral anomaly 

In the P-frame of the honeycomb, the Dirac equation of the two component fermionic 
wave \1/ = ("0, x) in the background field B is given by the following 2x2 matrix equation 




where T°, T^, are 2x2 gamma matrices which are related to the usual gamma matrices 
7^ as in (|2T|). 



1. Solving Dirac equation on supercell 

The gauge covariant derivatives Vi = di — iAi with Ai taken the SV gauge are as follows 

^^i = f >^^2 = ^ + ^Sx, (45) 

with no component Ax for the gauge vector. Using the expression of 7*^ = cr'^, 7^ = ia^, 
7^ = ia^ with the Pauli matrices as in eqQ, the matrices T* read as: 



where we have set ^ = uj\ — iu\ , ( = uf — iuj^ with the useful relations 
By replacing the weight vectors by their expression ffT^ . reads as 



(47) 



(48) 



CC = I , ?C - a = 

Notice that like in the C-frame, the Dirac operator T*Pj in presence of the background field 
B can be put as well into the form 

T^A = z I , (49) 



D 
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with D = + (1)2 and D = ST>i + C^^2- Its square, which is needed for solving the Dirac 
equation, is diagonal 



{VV,f = - I I , (50) 



and has the same eigen wave functions as T*Pj. Notice also that the explicit expression of 
D and D are as 



DD 



(51) 



and satisfy the commutation relation 

[AD] = -zi?(ec-ec) = -^, (52) 

where the sign of the right hand side depends on the sign of B] it is positive if sign(i?) is 
negative and vice versa. 



2. The explicit spectrum 
Setting 



D = J?1^A- , D = v/^At (53) 



the commutation relation fl52|) becomes 

[A-,A^ = , A^A- = A^A^ + ^/ (54) 

with the remarkable dependence in capturing the sign of the background field B. By 
substituting back into ( 150|) . the squared Dirac operator reads as 



DD 
DD 




(55) 



To get the solution of the Dirac equation, we first use the periodicity condition along the 
y-axis to expand the fermionic wave \l/ {x, y) in Fourier series like 

^(x,j/) = ^e^^^^z(x) (56) 

with 

(x) = I ^' ^"M , i^z. (57) 
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Substituting this expansion back into the Dirac equation, we obtain 





(58) 



where now the action of the operators D and D on the waves is restricted to the Fourier 
modes (x); this leads to 

Moreover, using the quantization property of the background field in terms of the area of 
the supercell namely BL1L2 = 2txQ, we can rewrite the above operators as 



(59) 



(60) 



Case 



-1 



In this case the operators 



-,D and A 



2\B\V^- — V 21^^ ^^"^ Heisenberg algebra 
El]) are respectively the creation operator and the annihilation one. So, we have 



(61) 



DD \ _2\B\V^ I A^A+l 
DD j 3 I A^A 

Thus, the wave functions \E'i = {ipi, xi) solving the Dirac equation are given by 



^/.n (x) 



4>l,n-l (x) 
<i)l,n {x) 



for the integer n > 1; and 



We also have, 



*z,o {x) 




0/,n {x) = ^ (A^y (pi^o (x) 

nl 



with the fundamental state (pi^ (x) obeying the condition 



d 



^— + <S (X+-L 



' dx 



I 



Q 



0i,o (x) = 



whose solution is given by 



i>i,o [X) 



oe 



2« 



(62) 



(63) 



(64) 



(65) 



(66) 



with C^=l-z^ and = |. We also have |0i,o| ~ e~l^l^('^+^^^)' 
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Case ]f| > 



In this case, the creation and annihilation are no longer as above; the new creation and 
annihilation operators obeying [C~ , = I are 



= J^,D , C- = J-i^D , (67) 



SO that 

f DD \ _,_^( C^C- ^ ^^g^ 



DD } \ etc- + 1 

]B\ 



The wave functions \E'; = (-0;, x'l) solving the Dirac equation with > are given by 



for the integer n > 1; and 



Knix)=\ I , (69) 



The excited states are given by 

<^'l,nix) = ^{CT^'lA^) (71) 

with ground state 

0;jx)=Aroe-^(^+*)(^-+^^O'. (72) 



3. Computing the index of Dirac operator 

Here, we focus on the case < and compute the index of the Dirac operator of fermions 
on a honeycomb supercell in presence of the background field B. Similar calculation can be 
done for the case p| > 0. 

• Case |f| < : Qtop < 

From eqs (l62tl63|) . we learn that the ground state $o (^j v) describing the fermionic wave 
with zero energy is antichiral 

^o(x,?/)= I ° ,a3^o = -^o. (73) 
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The function (f)g{x,y) is given by the following linear combination 

00 i^, V) = J]Qe^^^e-''(^+^^0' (74) 

where we have sei rj — \B\~ and where, a priori, the coefficients q are arbitrary moduli. 
As a wave function on the supercell, the function (f)^ (x, y) has to obey the periodic boundary 
conditions 

00 (x, y + L2) = (x, y) 
00 {x + Li,y) = e ^2 00 



The first condition giving the periodicity along the y-axis is manifestly satisfied due to 
f-e^iy+L,) ^ ^2^.l/f^^y^ The second 

requires the following identifications 



g*L2(^+^2) _ g2t7r«g* ^. The second condition giving the periodicity along the x-direction 



Q = Q_Q , / e Z (76) 

leaving afterward \Q\ free moduli 

ci, C2, c\Q\ . (77) 

This feature shows that the ground state 0^ (x, y) can be spanned as 

IQI 

•^0 y) = XI ^^^'"^^^0.- ^) ( 

with 

^o.(^,?/)-E^^^'^"'^^^*^^^^^' (79) 
Therefore the degree of degeneracy of the fermionic waves is as follows 



ground state : = , = |Q| 

excited states : iV+ = |Q| , = |Q| 



(80) 



So the number N^^^^ of degeneracy of the excited states is twice the number °^ of 
degeneracy of the ground state. The index of the Dirac operator which is given by — iV_ 
reduces to A^|'=°^ - 7Vl'=°^ = -\Q\ = Qt^. 

• anomalous QHE 
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To conclude this section, notice that applying the above results to the 2 Dirac points 
of graphene, we recover the well known value of the anomalous quantum Hall effect in 

2 

continuum. There, the transverse conductivity is given by oxy = ^aqhe\ with the filling 
factor v ^QHEi defined as the number of electrons per flux quantum, reads for the case of n 
Landau levels filled, as follows 

^.Q.. = 2(2n + l) = 4(^n+0. (81) 

Notice that the ground state corresponding to n = is half filled. From our analysis, we 
learn moreover that in the case < this ground state is filled by the zero modes \l/o 
eqs fl63ll73p having negative chirality; they obey the chirality property a"^\l/o = — ^^o with 
cr^ = diag (1, —1). In the case -j^ > 0, the ground state is filled by the zero modes \&o eq (l70|) 
with positive chirality satisfying cr^\E'Q = +^I^o. 



IV. FERMIONS ON HYPERDIAMOND 



In this section, we describe the supercell compactification in 4D hyperdiamond and study 
the solutions of the boundary conditions on fields required by periodicity properties. These 
tools will be used later on to derive the spectrum of lattice fermions on hyperdiamond in 
presence of uniform background fields B and E and to compute the topological index of the 
Dirac operator. 



A. Hyperdiamond and S\J (5) weight lattice 

Like 2D honeycomb, the 4-dimensional hyperdiamond can be described either by using 
the perpendicular compactification or the primitive compactification of the supercell. As 
the first compactification is a standard method in (hyper) cubic lattice QFT, let us focus 
below on the primitive one. 

• Hyperdiamond in P-frame 



The hyperdiamond L4 is a 4-dimensional extension of the 2D honeycomb; it is the world 
of fermions (light quarks) in lattice QCD4 j21|- 28|; see also 45|. In the primitive compact- 
ification, the hyperdiamond is generated by the basis vectors 



ai 



2 ^4 



(82) 
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= 2ui 


- W2 




012 


= 2UJ2 


- Ui 


- 0J2 




= 20J3 


- U2 


— U4 




= 2(^4 


-Us 





where a = dy | with d standing for the distance between sites and where cti, a2, as, are 
the 4 simple roots of SU (5). The reciprocal space L4 of the hyperdiamond is generated by 

bi = b, = '-^U2, b3 = ^^3, b4 = (83) 

where ui, U2, UJ3, oja are the 4 fundamental weight vectors of SU (5). We also have 



(84) 



• Unit cell of L4 

A unit cell in the hyperdiamond is given by the fundamental region made by the 4 basis 
vectors. In P-frame the basis generators are given by the vectors aj and the hyper-volume 
reads as 

ai A a2 A aa A a4 = — det (af ) (85) 

Like in the case of the 2D honeycomb, a unit cell contains two particle sites: one site inside 
the unit cell {the red ball in fig Hj) and 5 first nearest neighbors {sites in blue in fig H]) 
contributing each by a fraction | so that the total number isl + 5x| = 2. In the figure 
HI we give the projection of a 4D hyperdiamond unit cell on the 3D sub-space generated 
by {ai,a2,a3}; see also fig [3] for comparison with graphene. The first Brillouin zone of 
hyperdiamond is given by the unit cell in the momentum space generated by the bj's; its 
hyper-volume is 

647r^ 

bi A b2 A b3 A b4 = — — det (a;^) , (86) 

and is equal to with Vuc the volume of the real cell given by (185|) . Notice that the 4 
vectors bj (or up to a scale factor a;*) that generate the reciprocal space are dual to the a^ 
vectors (simple roots a^); they obey 

bj.aj = 27r6ij , a^.o;^ = afa;^ = 6^ . (87) 

and give a nice group theoretical interpretation of the relation between the real and momen- 
tum spaces. 

• More on hyperdiamond and SU (5) 
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FIG. 4: 3-dimensional projection of the real unit cell of hyper diamond in the primitive frame. It 
is generated by the 4 vectors ai , a2 , a3 , ( or up to scale factor by the 4 simple roots ai , Q!2 , 03 , 

04 of SU (5)). Sites Ai (resp Bi) belong to sublattice A (resp B) of the hyperdiamond. Bq has 5 
nearest neighbors of A-type. 

Hyperdiamond L4 is a 4- dimensional lattice contained in M^; it can be thought of as 
made of the superposition of two 4-dimensional sublattices A and B. Each site r„ in A has 

5 nearest neighbors (r„ + dei) belonging to B, 

r„ + (iei 

r„ + d&2 
r„ G A ^ ^ r„ + (/ea G B 
r„ + (ie4 
r„ + (ies 

and 20 second nearest ones r„ + (ej — e^) belonging to A. The 5 vectors e^, giving the 
relative positions of the neighbors, obey the traceless property 



ei + 62 + 63 + 64 + 65 = 0, 



(88) 



capturing the hidden SV (5) symmetry of the hyperdiamond. To parameterize the site 
positions r„ of L4, it is useful to use ei, e2, 63, 64 to generate vectors in hyperdiamond. 
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FIG. 5: On left the 5 first nearest neighbors in the 4^ hyperdiamond with the property ei + e2 + 
63 + 64 + 65 = 0. On right, the 3 first nearest in 2D graphene with ei + e2 + 63 = 0. 



These 4 basis vectors ej have the following intersection matrix 

/ 1 -i \ 



and can be realized as 



.1 1 _i 
4 4 " 








4 1 
-\ 



VI Vb VI I 

4 ' 4 ' 4 ' 4 



V^ -Va 



-VI 1 
4 



45414 

-v^ -VI V5 1 
4 ' 4 ' 4 ' 4 

-V^ Va -V5 1 

4 ' 4 ' 4 ' 4 



(0,0,0,-1) 



(89) 



(90) 



These 5 vectors ei, 62, 63, 64, 65 define the first nearest neighbors to (0,0,0,0) as in figO 
they satisfy the relations ej.ej = 1 and ej.ej = — ^ V i 7^ j showing that they are distributed 
in a symmetric way as required by the hidden SU(5) symmetry of L4. The 20 second nearest 
neighbors are given by 

and are remarkably generated by the 4 basis vectors 



a, = ¥ (e, - e2) , «2 = ¥ (^2 - es) ^^^^ 
03 = ^ (e3 - 64) , 04 = ^ (e4 - eg) 
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By using (pOl) . we have 



ai = (0,1,1,0) , as = (1,0,-1,0) 



as = (0, -1, 1,0) , ai= 



-1 1 -1 VE 

2 ' 2' 2 ' 2 



(93) 



from which we learn that af = 2 and det (ai, 02, a^, a^) = ^/5. We also have the intersection 



Ci' 



(94) 



matrix Cij = at.aj 

/ 2 -1 \ 

-12-10 
0-12-1 
\ 0-12/ 

indicating that the aj's are nothing but the 4 simple roots of SU (5) and C^j the correspond- 
ing Cartan matrix with inverse G*-' as 

/ 4 3 2 1 \ 

5555 



\ 



3 6 4 2 

5 5 5 5 

2 4 6 3 

5 5 5 5 

12 3 4 

5 5 5 5 



(95) 



This matrix G*-' is just the intersection matrix of the fundamental weight vectors oj* of SU(5) 



(96) 



with 



I 4||2 4 II 1||2 
kj = i , \\UJ \\ 
I II 5 ' II II 



\U! 



(97) 



1 i 1 

2' 2' 2' 10 

1^1 3VE\ 4 
2' 2 ' 2' 10 J ' 



The a;*'s span the momentum space of the hyperdiamond; their components in the canonical 
basis are 

(0,0,0,^^ 

with det (cj^, cj^, a;^, a;^) = — |v^- Therefore, generic positions r of and generic wave 
vectors k of the 4-dimensional momentum space decompose as follows. 



(98) 



X = u .r 



i=l 
4 

k = ^^fcjo;* , ki = aj.k 



(99) 



i=l 
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With these tools, we can determine several features on the crystallographic structure of the 
hyperdiamond; for instance the volume Vuc of the real unit cell is Vuc = ^V^', and the 
volume of first Brillouin zone Vbz reads as Vbz = ^5"V^- Recall also that the SU(5) Lie 
algebra has 10 positive roots and 10 negative ones (opposites); the positive ones consists of 
the 4 above simple and the 6 following 

a2 + as, a2 + a3 + a^, (100) 
cts + a4, ai + a2 + as + a^ 

Notice moreover that being also 4-dimensional vectors, the a;*'s may be decomposed as well 
like a sum over simple roots 

uji = |ai + |a2 + §"3 + 

5 5 5 5 ^^Q^^ 

^3 = |ai + |a2 + §"3 + |a4 
U4 = + |a2 + §"3 + |«4 

offering a unified way to parameterize the hyperdiamond L4; its sublattices A and B and 
their reciprocal spaces. 

B. Supercell compactification 

First we study supercell Csupcrceii compactification in hyperdiamond; then we move to 
deal with the boundary conditions of fields on supercell Cguperceii = SC4. 

1. p-cycles in supercell and 4-torus homology 

The supercell SC4 in hyperdiamond is 4-dimensional extensions of fig [3l and is described 
by a hyper-parallelogram V4 with periodic boundary conditions. Positions in supercell are 
parameterized by the vectors r = ^ . with the real coordinates = {x,y,z,T) com- 

pactified as 

< x' <li , li = aNi , x' = x' + li . (102) 

The 4D-parallelogram V4 is a real compact 4-cycle with the homology class of a 4-torus T"*; 
it has the following features: 
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has 4 basis divisors Aj and 4 dual 1-cycles Ei associated with the hyper-planes 
= in 7^4. These are: 





3- cycles 


duals = 


= 1-cycles 


X — 





A,. 


A, = 


Ex 


y = 





Ay 


A,= 


Ey 


z — 





A, 


A,= 




T — 





A. 




Er 



where Aj have the homology of a 3-torus and E^ of a circle S^. Positions in Ei and 
in Ai are given by the vectors 



El : r 
Ai :r 



X 



2 «1 

y ^OL2 + z r ^q;4 



(103) 



with the corresponding volume forms: 

Ex : dx ^OL\ 

Ai : ^dydzdr q;2 A cua A 

Similar relations can be written down for the others E^s and A^'s. 

has also six 2-cycles and six dual ones Cij given by 2D-parallelograms. 



(104) 





2-cycles 


dual 2-cycles 


X — 


y 


= 




— c 

^xy ^ZT 


X — 


z 


= 




Cxz — ^yr 


X — 


T 


= 


CxT 


CxT — ^yz 


y = 


z 


= 




c — r 

^yz 


y = 


T 


= 


r 


c — c 

^yr '^xz 


z = 


T 


= 


CzT 


CzT Cxy 



(105) 



These 2-cycles have the homology of 2-torii T^; they will be used later on to compute 
the fluxes of the background fleld components going through Cij. These 2-cycles are 
needed for the derivation of the topological index of the Dirac operator. Positions in 
these cycles are as 



Ci2 ■ r = X ^cti + y 
Ci2 -.r^ z ^a-i r 



(106) 
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with volume forms like 

Ci2 : \dxdy ai A 02 
Ai : ^dzdr a2 A 03 A 04 
Similar relations can be written down for the others. 



(107) 



the 4D-parallelogram V4, has 2^ vertices, they read in the P-frame as follows 

ihei, l2e2, he^M^A) , Q = 0, 1 . (108) 

with the liS as in ( 1T02D . 



Notice that as a 4-cycle, the 4D-parallelogram and its p-cycles are generated by the 1- 
cycle basis Ei, E2, E3, E4. Notice also that hyperdiamond sites in SC^ are given by discrete 
position vectors belonging to two classes 



respectively associated with the sublattice A and the sublattice B. Sites in A are given by 
the integer vectors 



= niai + n2a2 + n^a^ + n4a4 , < < A^j - 1 



(109) 



with SLi as in ( l82l) and n = (ni, n2, ns, n4). The sites of the sublattice B are given by the 
globally shifted ones, 



i-n = {ni + I) ai + (712 + I) a2 + (ris + f) as + {rn + |) 

with shift vector s = rj^ — given by 

s = |ai + |a2 + |a3 + |a4 
= d^y 

By substituting the vectors a^ by their expressions fl93|) . we obtain 

/o\ 






(110) 



;iii) 



:ii2) 



Notice also that the total number of sites in the supercell (I102p is 2N1N2N3N4] half of them 
belong to the sublattice A and the other half to B. 
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2. Boundary conditions 



To recover the hyper diamond lattice from a given supercell, one uses periodic boundary 
conditions on its borders as 

r = r + A''iai ,r = r + A^^^aq 

' ' ' (113) 

r = r + iV2a2 , r = r + iV4a4 

or equivalent ly 

r = r + k—ai. (114) 

To get the right spectrum of the Dirac operator, one has also to worry about the bound- 
ary conditions on the gauge field Afj, (r) and the fermionic waves (r) living on supercell. 
Generally these boundary conditions close up to a gauge transformation as follows 

Y(r + k^ai^ = (r) (115) 

where F stands for A^, and where fli some particular gauge transformation along the ctj- 
direction. To deal with the conditions (11151) . we extend the trick of Smit and Vink [5] that 
we have developed in section 3 by taking the gauge fields as 

AAr) = ,Ay{r) = -l3x ^^^^^ 
A, (r) = , Ar (r) = -Sz 

with no dependence in the y nor r variables. This particular gauge choice allows amongst 
others the two following helpful features: 

(i) it permits to use gauge invariance Af (r) = Ai (r) + diX to ensure the periodicity of 
the components of the gauge potential. This leads to a strong constraint on the y- and r- 
dependence into the gauge parameter A which are no longer arbitrary since they have to 
take the form 

X{r) = Bhy + ShT + <f{x,z), Q{r) = e'^^''^ (117) 

with some integration function {x, z) that we drop out below as it is not of great importance 
for the following analysis. Moreover periodicity of the gauge group element on supercell 
which has to obey i7 + ^^a^ j = Q (r) leads to the quantization of the background fields 
as follows 

Blil2=2nQB ,QBeZ 

(118) 

ShU = 27iQe ,QEeZ 
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(ii) it permits also to put the boundary condition on the fermionic waves into two classes: 
the two trivial ones 

^{x,y + z,t) = ^ (x, y, T, z) 
^ (x, y,z,T + k) = ^ {x, y, T, z) 
easy to deal with; and the two less trivial ones 



^{x + h,y,z,r) = {x,y,T,z) 
^{x,y,z + h,T) = {x,y,T,z) 

where fli and ^3 are particular gauge transformations to be given later on. 
Eq flll9l) allows to expand \1/ {x, y, z, r) in Fourier modes as follows 

,• 2irs 



(119) 



(120) 



^ {x, y, z, t) = J2e'~'J2^'~^ ^r,s {x, z) (121) 



restricting the problem of computing wave functions to looking for the Fourier modes 
^r,s {x-, z) whose dependence in x and z variables of ^r,s {x, z) can be obtained by solving 
the Dirac equation that we study below. 



V. SPECTRUM OF DIRAC OPERATOR 

Here, we compute the spectrum of the fermions on hyperdiamond. First we consider the 
case of the continuum limit (unit cell) . This is interesting to compare with the case of lattice 
fermions using supercell compactification to be considered later on. 



A. Dirac equation coupled to F^u- continuum case 



The Dirac equation describing the dynamics of four component Dirac fermions \E'^ (x), 
with SUc (3) color symmetry in presence of background vector potentials (x), is given 
by, 

443 

.r , 1 . . ^^^^^ 



where 



M=l /3=1 b=l 



a 
a 



ab 



mi = e^^ 

fi a 



a 
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and 



and where = ^^^""^^^ + Af^'^^^ with 

8 

^u.,Ai) = A;"^Q,m, = J^^r, (123) 

7=1 

with the 3x3 matrices 7} standing for the generators of SUc (3). Below, we focus on the 
particular case where 

^l^-(i) ^ , ^^^^('^ = , (124) 

but keep in mind that such analysis can be straightforwardly extended to include the other 
abelian components A^^"^^^ ^ leading to a non zero uniform field strength J^^u"^^^ along 
the Cartan directions in the SUc (3) Lie algebra. Notice that in (11221) the number e captures 
deformations of the energy spectrum around the Dirac point where live zero modes of QCD4 
fermions on lattice. Notice also that here we have ignored the flavor symmetry which is 
associated with the zero modes in the tight binding description of QCD on lattice; for 
details on implementation of flavor symmetry in the case minimally doubled fermions see 
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4^. 



• Dirac equation in continuum 

By setting Af,,^"^^'^ = in eq (ll22p . the Dirac equation in the canonical frame of the 
euclidian space reduces to its usual form, 

27^^ {d^ - lA^) ^ = e^. (125) 

In this equation, the 4x4 matrices 7^^ are euclidian Dirac matrices obeying the usual 4D 
Clifford algebra 7^'7'^+7'"7'^ = 26^'^ I4 with I4 the 4x4 identity matrix. These 7^ matrices can 
be built by using the Pauli ones r* and a* of the group product SUl (2) x SUr (2) ~ SO (4), 

with g"^ standing for both r* and a\ The g^, g^ satisfy the 2D Clifford algebra g^g' + g' g'' = 
25*-' J2 and the SU (2) symmetry bracket [g^, g"^] = 2ig^. We have 

(127) 

75 = ® , 7° = ® 
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with r^, cr^ = h 



I and 7^ = 7^7^7^7'*. More explicitly 
-ia^ \ .10 1 



7 



and 



r 



ia" 

/ 1 \ 

10 
0-10 

yo -1 y 



7 



/ 




, (g) cr^ 



/l o\ 
0-100 
0-10 
yo 1/ 



128) 



(129) 



— 7*7^ 



(130) 



The commutators j^^"^ give precisely the 6 generators of the spinorial representation of the 
SO (4) symmetry, 

where e*-'^ is the completely antisymmetric 3D Levi-Civita tensor. 

The Dirac operator (11251) involves also the gauge potential defined up to the gauge 
symmetry transformations 



with A (x) arbitrary gauge parameter. Since the background field strength F^^ 
di,A^ is a constant, the gauge potential reads as follows. 



(131) 



Notice that generally F^i, has six real degrees of freedom, 3 magnetic Bi and 3 electric Ei. 



flU 



( +^3 -S2 

-S3 +5i +^2 

+S2 -Si +S3 

-^2 -^3 / 



(132) 



• Choice of the background fields 

Below we will think about this F^^^, as a sub-matrix of the following 5x5 antisymmetric 



one, 



F F K 




(133) 
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where now the F^s's are the 4 components of the electric field in 5D; and F^j, are the 6 
components of the magnetic tensor. We also make the two following useful choices: 
(i) we restrict the matrix to the particular case, 

/ \ 

f5 
-F 
+E J 

allowing exact computations due to the decoupling of the left and right sectors of fermions. 



\ 



(134) 



Using the antisymmetric tensor Snupa, we have 



(ii) To deal with the gauge field, we can either use the symmetric choice 



9 ^ ) 



A, 



A. — -^7 



(135) 



(136) 



or the Smit-Vink method 



(137) 



Ai = 0, A2 = Bx 
A3 = 0, A4 = Ez 

In the next sub-subsection, we use the first choice as it allows to take advantage of the 
symmetric role played by the components fields to solve the Dirac equation in the continuum. 
Later on, we use the second choice to study the Dirac operator of fermion on supercell; the 
gauge f ll37p is convenient for the study the degeneracy of the zero modes of the Dirac 
operator and its topological index. 



B. Spectrum in the gauge ()136p 



To get the spectrum of the Dirac operator ( 1125^ . notice that the 4 gauge covariant 
derivatives Di, D2, -D3, -D4 satisfy the generic commutation relations 



but because of the choice (I134p of the background fields, they reduce to, 

[DuD^] =iB , 



[D3,A] =iE , 



(138) 



(139) 



and all others vanishing. 
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1. Deriving the wave functions 



The fermionic wave functions solving the Dirac equation are given by representations of 
the algebra f ll39p . By using f ll36p . the four covariant derivatives organize in 2 + 2 capturing 
a complex structure as follows, 



D,-^D, = f^ + ^u,D, + ^D, = f^-9iu 



(140) 



In getting these relations, we have used the explicit expressions Di = di — 'i^Ui -^2 
82 + together with similar relations for D3, D^; and we have set 



The representations of eqs fll39p depends on the sign of B and E. Setting 

i (Di - iD2) = A-y^2\B\ ,i{Di + 1D2) = A+y^2\B\ , 
I (D3 - iD^) = C-./2\E\ ,t{D3 + iD^) = C+./2\E\ , 

the commutation relations f ll39p read also as 

[A , A+] = , [C , C+] = |f| / , 
[A-,C^]=0 ,[A+,C±]=0 ; 



(141) 



(142) 



(143) 



These relations f ll43p show that the Dirac fermion in the background field F^^, fll34p describe 
a priori 2 quantum harmonic oscillators with oscillation frequencies 



^/2\B\ , z/' = ^/2\B\. (144) 



The operators A'^A and C^C , which a priori give the number of energy excitations in u 
and z/' units respectively, read in terms of the gauge covariant derivatives as follows 



2 \B\ A+A- = {D,y + {D2r + t[D,, D2] 
= {D,f + {D2f-B 



and similarly 



(146) 
{D^f + {D,f - E 
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Using the expression of the matrices 7^^, we can write this matrix operator Ha/3 — ^la^Df^ 
in terms of the "creation" operators C+ and the "annihilation" ones C~ as follows: 

/ -i/C+ -uA-\ 

10 -vA+ u'C- 

i v'C- vA- 

^ -i/'C+ y 



H. 



a/3 



(147) 



with — H. Moreover, using the commutation properties C'^A^ — A^C'^ and C — 
A^C", the squared hamiltonian H'^ reads as follows, 

/ u^A-A++u'^C+C- \ 

iy^A+A-+iy'^C'C+ 

u^A-A++i/^C~C+ 

^ iy^A+A-+iy'^C^C- J 



and, by using = {ipa,id), leads to 



with 



00^ 



^2 / \ ^2 



6 / I 6 



ly^A+A- + iy'^C-C+ 

u^A-A+ + u'^C-C+ 

u^A+A- + i/'2c+C- 



(148) 



(149) 
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showing that the two Weyl spinors ipa and can be treated separately. Putting the expres- 
sion of CO^ and O^O back into eq ffTiHj) . we obtain 



(u'^C-C+ + u^A-A+)^, =e^^i , 



(150) 



The solutions of these equations depend on the sign of the background fields B and E. In 
the case > 0, -j^ > 0, we have: 



V6/ 



On, m—1 

©n— 1, m—1 
Or),_ m. 



(151) 



with 



and 



0n,m {u, U, V, v) = 6'„ {u, u) X 6'^ (f , ^;) , 

A-xeoiu,u) = {§^ + ^u) 00 iu,u) 



ra, m > 0, 



(152) 



(153) 



solved as follows 



v) e i 



(154) 



where the complex functions Aq (w) and Mq {v) are anti-holomorphic functions. The excited 
waves functions On {u, u) and O'^^ (f , v) are obtained by applying the creation operators. 



(155) 



2. Zero modes and topological index 



The zero modes of the Dirac operator depend on the sign of the background fields B 
and E since the algebra of the commutation relations f ll43p depend on and as given 



\E\ 



below: 



[A , A^] — , [C , C+] — . 



\E\- 



(156) 
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We have the 4 following possibihties: 

(a) Case ^ > 0, ^ > 
In this situation, the A'^, (7+ are creation operators and A', C~ annihilation ones. So the 
zero mode is given by 



e=0 



/ \ 




V 00,0 J 



(157) 



with the chirality property 



(158) 



(0-3 O T^) *,=o = +*£=0 
7^*6=0 = -*e=0 

and by setting Tr (^£=o^£=o) = 1, we have 

Tr [f'e=o7'*e=o] = -1 

(b) Case Ji < 0, ^ < 
In this case A'^, C"*" are annihilation operators and A~ , C~ creation ones so that the zero 
mode is 

/ ©0,0 \ 



(159) 



*6=0 






V / 



(160) 



and 



as well as 



{a^ (g) r^) ^,=0 = +*£=o 
7^*6=0 = +*£=o 

Tr [^,=Q {<j^ ® T^) *,=o] = +1 
Tr [^,=o7'*e=o] = +1 



(161) 



(162) 



(c) Case ^ > 0, li < 
Here , C~ are creation operators and A~, C~ annihilations. The zero mode is given by 

/ \ 



©0,0 

V > 



e=0 



(163) 
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with 



and 



7^*6=0 = -*e=0 



Tr [*,=o7'*e=o] = 1 

(d) Case ^ < 0, ^ > 
In this case, the zero mode reads as 

/ \ 



e=0 



with 



and 



©0,0 



V ) 



(cr^ (g) r^) *,=o =-*£=o 
7^*6=0 = +*£=o 

Tr [t,=o (8) T^) ^^,=0] = +1 
Tr [^,=o7'*e=o] = -1 



(164) 



(165) 



(166) 



(167) 



(168) 



Prom the above analysis on the chirahty of zero modes, it follows that the chirality operator 
that satisfies the Atiyah- Singer theorem is ® since 



Tr [^,=0 {o^ ® T^) *.=o] = 1 



(169) 



This feature can be explained as due to the factorization of the SO (4) symmetry of IR^ 
in terms of the product of S'f/i (2) x SUr{2). The relation that involves 7^ reads as 
Tr [^',=o7'*.=o] = ^. 



VI. SOLVING DIRAC EQUATION ON SUPERCELL 

To study the spectrum of the Dirac equation in the supercell compactification of the 
hyperdiamond, we extend the result of sub-section 5.2 concerning the continuum to the case 
of the 4D lattice. There, we have used the canonical frame of with the local coordinates 
— {X, Y, Z, T); which will be used also in the case of lattice. 
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In the C-frame, the 4 gauge covariant derivatives = {Di, D2, D3, D4^) satisfy the usual 
commutation relations giving the components of the U(l) gauge curvature 

[D^,D,] = -iF^,. (170) 

By choosing the background fields as in (11341) . these commutation relations factorize into 
two decoupled Heisenberg algebras as follows, 

[Di, D2] = iB , [Di, D,] = [Di, D,] = ^^^^^ 
[D3, D,] = lE , [D2, D3] = D^] = 

To work out explicit solutions of these equations, one may used either the potential vector 
or (11371) . In this section, we use eq (ll37p leading to the following gauge covariant 



derivatives = — iA^ 

Di = , Do = — iBX 

■^3 = T D4 = Qj; — lEZ 

This choice breaks the symmetric role played by the variables {X,Y) and {Z,T); but is 
suitable to deal with the boundary conditions of the fields on supercell. 
Notice that eqs (ll70l) and (I172p can be also expressed in the P- frame with positions as 
X* = (x, y, z, t). The passage between C- and P- frames is given by the transformations 

= afx^ , x' = u^Xf^ , a':ui = 5l (173) 

Similar relations can be written down for the other the objects; for instance the potential 
vector Afj^ and the gauge covariant derivatives in C-frame are related to their homologue 
Ai and Vi in the P-frame as 

^M = ^M- D, = ujlP, ^^^^^ 



A. Computing the fluxes of background flelds 

The flux of the background fields through hyperdiamond supercell is a scalar quantity 
and is frame independent. This flux give the total topological charge inside the supercell 
SC4; it controls the chirality of the ground state and allows to determine the topological 
index of the Dirac operator in the background fields B and E. 
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To compute the flux, one can either use the C-frame or the P-frame; in fact it is frame 
independent. To see this property recall that in the C-frame the gauge curvature is given 
by F^i, and in the P-frame is J^ij: 

[D^, D,] = -iF^, , [Vi, Vj] = -iTij . (175) 
These two tensors are related as 

J u.^ u.^± 2 «j ; fJ.f (176) 

with 

ai A ttj = (af - a^^a\) 

, . . (177) 

The corresponding gauge invariant 2-form fleld strengths are then given by 

^ ^ 178 

F = IdX^' A dX'^F^^ 

and are equal = F since they are frame independent; thanks to afw^ = Sf. Moreover, 
because of the choice 

F^u = -B£^iy34 + Eeuflu (179) 

that we have used in this paper to work out explicit solutions of the Dirac equation, the 
2-form F reduces to the simple expression 

F = BdX ^dY + EdZ A dT (180) 
By using the coordinate change fll73p to the P-frame, it can be also written as 

F = [Baja] + Eafa"^) dx' A dx^ (181) 

where F takes a general expression in the basis dx^ A dx^ . With these relations at hand, we 
can compute the flux of the background flelds through the various p-cycles of the supercell. 
We will do the calculations the C-frame. 

The total topological charge Qtot of the background fleld within the supercell is given by the 
integration of the 4-form F A F over the supercell, 

(^1^ /" ^FAF = Qtot , Qtot e Z , (182) 
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Substituting fllSOp back into fll82p . we obtain a quantization condition on the background 
fields given by 

BxE = j \ ' Qtot (183) 

One can also compute the fluxes 

tJj ''''' 

of the fleld strength F through the 2-cycles Cij = C2 of the supercell; this gives extra 
quantization conditions. Because of the choice (11341) . non trivial fluxes are indeed given by 
the 2-cycles C12 and C34. Using the relation F = dA to map the integration over the 2-cycles 
Cij to circulation around its boundaries dCij, we end with 

j_ / r _ 1 

JCu JdCi2 

and 



F = ^f A = Qb (185) 



I F = lf ^ = Qe (186) 

JC34 JdC-ji 

giving two extra quantization conditions; one on the background fleld B and the other on 
E. More precisely, we have: 

L1L2 



B = y^Qb ■, Qb 



(187) 



= T^fi^ 1 Qe 



Comparing (11871) with (I183p we get the following relation between the topological charges 

Qtot = QbQe- (188) 



B. Dirac equation on 4D- supercell 

The euclidian Dirac equation on supercell is given by the 4-dimensional extension of the 
Dirac equation on honeycomb lattice. In addition to periodic background potentials, this 
equation involves four component fermionic waves with boundary conditions described by 
the SO (4) spinor \Ef = {^a^^aj that we want to determine below. 



1. The hamiltonian 



The euclidian Dirac equation is 

HI - \=e \ - \ (189) 
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where the hamiltonian H reads in the canonical frame as H = 'y'^D^ or equivalently in the 
P-frame 

H = = i^J;d, = rVi (i90) 

with r* = The solution of f ll89p in continuum has been worked out in the previous 

section; below we want to extend these results to the lattice case where boundary conditions 
put strong constraints on the solutions. To that purpose, let us start by collecting useful 
tools. First the hamiltonian H has the form 



with 



D 



H 



\E\C- 




(191) 



and 



D 



and where we have set 



A- 



'\B\ 



-iJ\B\A' 
\E\C+ 



\E\C+ i^/\B\A- 
t^/\B\A+ ^/\E\C' 

[D, - iD^) , C- = 



(192) 



(193) 



(Di + iD^) , C+ 



(D4 + iD^) 



(194) 



/\B\ ^ ^ ^' ' ^\E\ 

Second the operators A"^ and obey the commutation relations 



(195) 



and all other vanishing. Notice that these commutation relations have a remarkable depen- 
dence on the sign of the background fields B and E. Third, the squared hamiltonian has 
the diagonal form 

^ DD 
DD 

with 

^ \B\A-A-^ ^\E\C-C+ 
^ \B\A+A- + \E\C+C- 



DD 



(196) 
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and 

I \B\A-A+ + \E\C+C- 
= ' ' ' ' (197) 

\ \B\A+A- + \E\C-C+ 

involving the four possible quadratic combinations of and namely A^v4+, A'^ A~ , 
C-C+ and C+C'. 



2. The solutions and Index{D) 



The solutions of (11891) on supercell are representations of the algebra (11951) that have to 
satisfy moreover the boundary conditions on the fields. These conditions are quite similar 
to those studied in the case of 2D honeycomb; so we omit here the lengthy technical details 
and just give the results. 

There are 4 classes of solutions of (I189P depending on the signs of and j^. They are 
obtained as follows: First, expand the wave ^ (x) on the periodic supercell in Fourier series 

as 

V * ■ f 2liT I 2q-7r \ 

^'^^^y^^.^)m,^^{x,z) (198) 



* (x) = Y^e' 

l,q 

This expansion follows from the periodicity of \E' (x) along the y- and r- axis. Second, solve 
the non trivial boundary conditions along the x- and z- axes by following the method used 
in the case of 2D honeycomb which lead to eqs (l42tl43|) . As in the present case we have to 
deal with the 2 variables x and z, we write the {x, z) modes like 



(199) 



with 



6 



X 



z + 



Qe 



(200) 



The next step is to determine the function $ (^, this is a Dirac spinor which we set as 

1 0(?,C)\ 



(201) 
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(202) 



with components obeying the following coupled equations 

{\B\ A-A+ + 1^1 C-C+) (j) = e^(j) 
ilBlA+A- + \E\C+C-)g = e^g 
{\B\ A-A+ + \E\ C+C-) if = eV 
i\B\A+A- + \E\C-C+)x = e\ 

The operators A"^ and satisfy the algebra (I195p and show that solutions for $ () 
depend on the sign of the background fields. These solutions are as follows: 

• case > 0, ^ > 

In this case, which corresponds also to a positive topological charge Qtot, the algebra 
f ll95p reads as follows 

\A+,A-] = 1 , 

^ ' ^ ' (203) 

showing A^, C~ are creation operators and A^ , are annihilation ones. Using general 
results on quantum harmonic oscillators and the relations 

Ma = A-A+ , A+A- = A-A+ + 1 
Uc = C'C+ , c+c- = C-C+ + 1 

where M stands for the number operator, it is not difficult to see that the e energies are 
discrete as 



=n\B\+m\E\, n.mEN 
and the corresponding wave functions are given by 

/e \ 



(204) 



e 



n — l.m — 1 



0. 
e 



(205) 



where 



and 



en,miC,O = 0n{Ox^n. (C) 

9n{0 =^{A-reoiO ,A+9oiO = 

^mio = ^ {c-r ^oio ,c-^, (0 = 



(206) 
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with 



^0 (0 = A/'oe-^«\ ^0 (C) = AToe- W 



(207) 



Notice that because O^i = 'd-i {() = 0, the ground state has only one component as 



'0,0 





(208) 





V J 

Moreover since the degeneracy of 9o and do are respectively \Qb\ and IQ^I; it follows that 
the degree of degeneracy of $o,o is precisely the total topological charge 



\QbQe\ = \Q: 



tot\ 



(209) 



in agreement with Atiyah-Singer theorem in 4-dimensional Dirac theory. 



case > 0, ^ < 



This case corresponds to Qtot < 0; the algebra (11951) reduces to 

[A+,A-] =1, 
[C',C+] = 1 . 



(210) 



It shows that A , C"*" are the creation operators and A'^, C are the annihilation ones. The 
energies ^ are same as above but the fermionic wave are like 



e 



\ ^n-l,m-l / 

Here also the ground state has one component given by 

/ \ 



©0,0 

V / 

it has the same degree of degeneracy as in the previous case. 



(211) 



^ n 



(212) 
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B 

\B\ 



case < 0, |f| > 0, with Qtot < 



(213) 



The algebra fll95p becomes 

[v4-,A+] =1 , 
[C+,C-] = 1 . 

Here A'^ , C~ are creation operators and A~, are annihilations. The fermionic waves are 
as follows: 



(f ■ C) 



1 1 

Ti — 1 , m — 1 



(214) 



with ground state as: 



/ \ 




V ©0.0 / 



(215) 



case 



B 



1^1 < 0, ^ < 0, Qtot > 
The commutation relations fll95p 

[A-,A^] =+1, 

[c-,c+] = +i, 

with A'^, the creations operators and A^ , the annihilators. The fermionic waves are 



(216) 



with ground state like 



0,0 



are as 




/e 

n- 
















^ 




©0,0 









V 



n — l,m — 1 



\ 



(217) 



(218) 



The index of the Dirac operator is given by Tr \^e=o {cr^ ® t^) ^e=o] = Qtot- 
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VII. CONCLUSION AND COMMENTS 



In this paper, we have studied topological aspects of fermions on a family of 2N- 
dimensional lattices in presence of background fields with special focus on the 2 leading 
crystals namely the graphene and the 4D hyperdiamond of QCD^. With the results ob- 
tained by our explicit study, we have now an exact answer on the population of the ground 
state of fermions on lattices in presence of uniform background fields. For example, in the 
case of graphene in a strong magnetic field, we find that the chiral anomaly is behind the 
observed anomalous in the filling factor Vgra = 4 (n + |) of the integer quantum Hall effect. 
This means that the ground state of graphene with Vg"^^ = 4 x ^ = 2 is occupied either 
by positive chiral states or negative ones depending on the sign of the magnetic field B. 
The same statement can be made for light quarks of QCD^ on hyperdiamond and more 
generally fermions on higher even- dimensional honeycombs. In QCD4 on lattice with a 
Dirac fermion (say the quark u) in the background fields B and E; the filling factor reads 
as lyqcD = # ('^ + I) ("^ + i) there are 4 possible configurations for the population of the 
ground state depending on the signs of B and E; in the case of u quark these are u^l, uil, 

UtR, UlR- 

To exhibit this behavior, recall that in honeycomb the magnetic field B appears in the 2 
gauge covariant derivatives Di and D2 whose curvature can, up on using the scaling ( I53ll54l) . 
be put into the remarkable form 




This is a typical Heisenberg algebra; but with two sectors depending on the sign of the 
magnetic field B. For positive B, the operators and A~ are respectively the creation 
operator and the annihilation one; but for negative B, this property gets reversed since now 
A~ plays the role of the creation operator and A'^ the role of the annihilation one. Almost 
the same thing happens for quarks in QCD4 and fermions on higher dimensional lattices. For 
example, in the case of fermions on 4D hyperdiamond, we have 4 gauge covariant derivatives 
Dfj^ = {Di, 02,0^, D4) obeying the general commutation relations 

where F^y is a C-number capturing in general 6 degrees of freedom (provided < F^^^ > 
= 0). A careful analysis of this algebra shows that it describes two interacting quantum 
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harmonic oscillators. However by choosing the gauge field strength as Ffj,^ = Bie^u^i + 
B2Si2fj.u, the above algebra reduces to 2 uncoupled Heisenberg algebras 

but with 2^ sectors according to the signs of Bi and B2. This result extends straightforwardly 
to the case of fermions on 2 N- dimensional honeycombs with background fields. There, the 
field strength F^i, has generally {2N — 1) moduli and so the corresponding algebra of 
the gauge covariant derivatives describes interacting quantum harmonic oscillators. By 
choosing the field strength as 

Fi_i,y = -Bi 6^1^34.... 2 AT + -82^12^41^56.. ..2Af + ••• + B]\f 612... ./j.^, 

the algebra of the covariant derivatives gets reduced to N uncoupled Heisenberg ones as 
above and has 2^ sectors depending on the signs of the Sj's. In this generic case, the total 
topological charge Qtot of the background fields within the compact supercell SC2N of the 
2 N- dimensional honeycomb as well as the partial fluxes through the 2-cycle basis Ci of the 
supercell are given by 



a 



Qtot = Kr / F A ... A F, Qi = — 

rN 



niV 
Qi which can be proved by thinking about the supercell 
i=l 

SC2N as given by the product of those 2-cycles Ci of SC2N with no intersection; i.e: CiHCj = 
0. The computation of these charges for the case N = 2 was done in section 6; see eqs fll87t - 
11881) : they can be easily extended to higher dimensions; in particular for the case fermions 
on the 6-dimensional honeycomb. By taking = BiS^y^Am + B2ei2^ub%+ B^ei234fj,u, one 
ends with Qtot = QiQ2Qz- 

In the end, we would like to add that our explicit approach gives as well a unified group 



theoretical description of fermions in both graphene and QCD4. The construction of 20 1 
turns out to be intimately related with the weight lattice of SU (3) we have given in section 
2 and 3; the primitive compactification used in [20] has also to do with the simple roots 
basis of the SU(3). The latter is a hidden symmetry of the 2D honeycomb; it allows many 
explicit calculations in graphene and moreover draws the path to follow to get the extension 
of results on graphene to fermions on 2 N- dimensional honeycombs where the job is done by 
the hidden SU {2N + 1) symmetry. 
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VIII. APPENDIX: LATTICE CALCULATIONS 



In this appendix, we give some details on the lattice calculations used in this study. These 
computations, which have been understood in the paper; are based on the method of refs 
% |5| ; and constitute an extension of results, obtained in ref 20| concerning graphene, to 
the case of QCD on 4D hyperdiamond. For completeness, we describe below the 2 following 
useful things: 

(1) review briefly the main lines of tight binding model for graphene; a quite similar analysis 
is valid for QCD4 on hyperdiamond (see section 4). We also take this opportunity to develop 
further the link between the electronic properties of graphene and SU (3) representations. 
This group theoretical approach extends directly to 4D lattice QCD on hyperdiamond. 
There, the role of SU (3) is played by SU (5). 

(2) give a comment on the reason behind solving Dirac equation with periodicity conditions 
on the wave functions. This is important for two things: (z) working out the wave functions 
with explicit dependence on lattice parameters, the quantized topological charge Q; and too 
particularly the exact determination of the degree gm of degeneracies of the states of energy 
Em as required by the computation of eq(IT]). (u) performing numerical calculations as done 



in 20| for 30 x 30 lattice to check analytic predictions on the Atiyah-Singer theorem, the 



chiral anomaly and related issues. 

A. Case of graphene: a brief review and link with SU (3) 

On the 2D honeycomb lattice, made by the superposition of two sublattices A and B as 
depicted by figs [T]|2l the tight binding hamiltonian describing the hopping of electrons of 
graphene to first nearest neighbors, in presence of a magnetic background field B, reads as 

H=-tJ2 ( Wr„,v. feX.+v. + W;„,v, <-vJ (219) 

rmGA \ 1=0 1=0 J 

Here t is the hop energy; af^ and &^„+v, are the fermionic waves respectively associated with 
A and B sublattices and satisfying the usual anticommutation relations. Wr„,vi is the link 
field given by 

^rm,v, = e^^'-^('"'") = e^'''"-^^ (220) 

with Af^ [x, y) the potential potential of the external magnetic field. This is a u (l)-valued 
gauge connection emerging from the site Tn and lying along the v/ direction. Notice that 
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if switching of the external B, the field Wr„,v, reduces to identity; and then H describes a 
tight binding model of hopping free electrons. We also have the following relations 

Vi.^ = +lAi + 

V2.^ = -H + (221) 

Vq.^ = -|A - f A 

where Ai stands for a^A^ (r) with cti , 02 the vectors generating the sublattices A and B of 
the honeycomb. In our present study, we have used the method of 0, 0] to deal with the 
gauge potential A^ (r) by working with the helpful gauge choice 

Ai = at A, (r) = , A2 = ai^A, (r) = -By (222) 

Recall that on honeycomb each pi-electron of a carbon atom, say a fermion ar„ of the 
sublattice A, has 3 first nearest atom neighbors of B-type with fermions &r„+vo) ^r„+vi, &r„+v2- 
The vectors vq, vi, V2 with entries v; = {v^) parameterize the relative positions of &r„+v; 
with respect to ar„ and satisfy some remarkable features; in particular the following: 
(a) they obey the vector constraint relation 

vo + vi + V2 = (223) 

that captures physical information on the phases of wave functions. 

(6) they play a central role in dealing with lattice calculations as they encode the hopping 
to nearest neighbors. In particular, they allow to build the analogue of the usual derivative 
term of the Dirac hamiltonian in continuum; the lattice derivative turns out to be captured 
by some combination of the phases e'^^'^'-. 

(c) in our approach, the v^'s are interpreted as the 3 weight vectors of the fundamental 
representation of the SU (3) group. This is a remarkable observation that allows to simplify 
drastically the lattice calculations. To get the point, set 

Vo = Ao , vi = Ai , V2 = d A2 (224) 

with d ~ 1A2A the length of the carbon-carbon bond in graphene; then put back into (12231) . 
we end with a well known group theoretical identity namely 

Ao + Ai + A2 = (225) 

describing precisely the constraint equation on the weight vectors of the SU (3) fundamental 
representation. 
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Recall by the way that for SU (3), the weight vectors of the complex 3-dimension represen- 
tation can be written in two basic manners: either in terms of the fundamental weights oji 
and UJ2 of SU (3) as follows 

Ai = +uji , A2 = -ui +UJ2 , \o = -0J2 (226) 
or equivalently in terms of the two simple roots ai and 02 of SU (3) like: 

Ai = +!«! + |a2 

A2 = + |«2 (227) 

Aq = ~ §Q?2 

satisfying manifestly f l225p . Obviously the fundamental weights wi, uj2 and the simple roots 
«!, a2 are linked by the duality relations ai.ujj = 6ij which lead to ai = 2uji — UJ2, 0.2 = 
2u2 — Wi- We also have the helpful relations 

«! = Ai — A2 

«2 = A2 — Ao (228) 
tti + tt2 = Ai — Ao 

in dealing with lattice calculations. This SU (3) group theoretical analysis has been shown 
to extend straightforwardly to the 4-dimensional hyperdiamond of section 4 where the role 
of SU (3) is played by SU (5); see eqs(l82D-(IT0Til. 

If switching off the interaction between the electrons of graphene with the external magnetic 
field B\ and then performing the Fourier transform of the local fields, we can put the tight 
binding hamiltonian H\s=q into the form Xlk-^k with wave vectors k = [kx-, ky) and 




H^.= [a-^■X) I _ I I I (229) 

with 5k = 'Ylii e^'^^'^'- or more explicitly 

£k = e''^^-^° + e''^^-^' + e'"^^-^^ . (230) 

Moreover, using eqsf l2261l228p . we can also put the above relation into the remarkable form 

which is manifestly invariant under translations in the reciprocal lattice; that is under the 
shifts 

k k+^(r2ia;i + n2a;2) , r2i,n2 G Z (232) 
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The diagonalization of the hamiltonian mode ifk fl229p leads to the energy dispersion rela- 
tions 

E± = ^'^'t\j ^ 2 '^^^ d\<i.ai + - cos (ik.a2 + - cos dk.cta (233) 
where we have set = ai + a2- 

Notice that zero modes of the hamiltonian Hy_, obtained by solving the vanishing condi- 
tion £k = 0, are immediately learnt from (123 ip and are given, modulo translations in the 
reciprocal lattice, by the two following Dirac points 

k± = ±^ (wi + U2) = ±kF (234) 

Notice moreover that near these zeros, say for k_|_ = q + k^? with small q, the dispersion 
energy relation is linear in q and the physics of the electron is mainly described by a free 2- 
dimensional Dirac theory in continuum with the periodicity property (12321) of the reciprocal 
lattice. 

By switching on the interaction between the electrons and B, the previous energy dispersion 
relation = (^d.x'^k^, g^^g jxiodified; it is given by a complicated expression which is 
obtained by substituting the wave vector A;^ by the gauge covariant quantity K,^ = k^—A^. 
Furthermore, using the fact that the external magnetic field is constant, we end, after some 
straightforward algebra, with the following covariant derivative in reciprocal space 

]C, = k,- ^^e,,^^ = A + _^-p/,^^ (235) 

leading in turns to e (K.) = ^^e^d.K'^i^, Near the zeros (I234p . e (K.) leads therefore to the 
Dirac equation in continuum in presence of a background field B. This equation reads in 
the real space as, 

- ^ix,y) = 0, ^ = I M (236) 




with the periodicity property 

r r + (niOi + n2«2) , ni, 72-2 G Z (237) 

In section 3, we have worked out the solutions \1/ (x, y) near the Dirac points that obey 
periodicity properties eq( l39|) and ( 129H3T]) . It was shown that the ground state is chiral and 
depends on the sign of p-. In the case -j^ = — 1 for instance, the wave functions are given 
by eqs fl^ . flMl) and where, in addition to the background field B; the dependence 
on the lattice parameters, the topological charge and the degeneracy of the eigenvalues are 
manifestly exhibited. 
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B. Numerical study 



In order to analyze the spectrum of the Dirac operator D of graphene in various gauge 
field backgrounds B = j^^, one diagonalizes the matrix D'^ by using a subspace iteration 
technique as well as Chebyshev polynomial iteration to accelerate the convergence of the 
eigenvalues. Following {20I; see also j^] for technical details, the plot the 60 smallest 
eigenvalues E"^ (1 < m < 60) reveals a very good agreement between the numerical results 
and the analytic predictions as shown by fig El These energies are calculated on a A'^i x N2 
lattice supercell (A^i = N2 = 30) with primitive boundary conditions and for values of 
topological charge Q varying between 1 and 26. 




Recall that the analytic prediction for the eigenvalues of the energy spectrum leads to 
= ^m\B\ with m a positive integer. Lattice calculations show also that the degeneracy 
pattern (139]) of these eigenvalues is as follows 

J \Q\ for m = 
gm= < (238) 
1 2 \Q\ for m > 

By substituting the background field B by its expression in terms of the topological charge 
Q and the area of the supercell, we can rewrite the energy spectrum like, 

E"m = ^^\Q\m (239) 
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As shown on the plot, the states with Eq = 0, that form a triangle on fig [6l corresponds 
precisely to the zero modes of the Dirac operator with degree of degeneracy gm growing 
linearly with Q in complete agreement with the index theorem. 

Acknowledgement 1 : 

The authors thank the Hassan II Academy of Science and Technology where part of this work 
has been done. E.H.S thanks the Moroccan Center for Scientific Research and Technology; 
Project ref URAC09, for support. 



[1] M. Atiyah, I.M. Singer, Ann. Math. 93 (1971) 139, 

[2] M. F. Atiyah, R. Bott, V. K. Patodi, On the heat equation and the index theorem, Inv. Math. 
19 279 (1973), 

[3] I. Barbour, M. Teper, Phys. Lett. 175B (1986) 445, ELSEVIER. DOI: 10.1016/0370- 
2693(86)90621-0, 

[4] J. Smit, J.C. Vink, Phys. Lett. 194B (1987) 433, DOI: 10.1016/0370-2693(87)91078-1, 
[5] J. Smit, J.C. Vink, Nucl. Phys. B286 (1987) 485, DOI: 10.1016/0550-3213(87)90451-2 
[6] S. Itoh, Y.Iwasaki, T. Yoshie, Phys. Rev. D36 (1987) 527, 
[7] J.C. Vink, Nucl. Phys. B307 (1988) 549, 
[8] T. Kalkreuter, Phys. Rev. D51 (1995) 1305, 

[9] W. Bardeen, A. Duncan, E. Eichten, H.Thacker, Phys.Rev. D59 (1999) 014507, 



hep-lat/9705002 ; Phys.Rev. D57 (1998) 1633-164, |hep-lat/9705008 , 
[10] David H. Adams, Index of a family of lattice Dirac operators and its relation to the non-abelian 

anomaly on the lattice, Phys.Rev.Lett.86:200-203,2001, |arXiv:hep-lat /9910036] 
[11] F. Karsch, E. Seller, I.O. Stamatescu, Phys.Lett. B157 (1985) 60 , DOI: 10.1016/0370 - 

2693(85)91212-2, 

[12] Keun- Young Kim, Bum-Hoon Lee, Hyun Seok Yang, Zero Modes and the Atiyah-Singer Index 



in Noncommutative Instantons, Phys.Rev. D66 (2002) 025034, |arXiv:hep-th/0205010 
[13] Ting-Wai Chiu, The Index and Axial Anomaly of a lattice Dirac operator, 

Nucl.Phys.Proc.Suppl.l06:715-717,2002, |arXiv:hep-lat/0110083t 
[14] David H. Adams, Families index theory for Overlap lattice Dirac operator. I, Nucl. Phys. B624 



(2002) 469-484, arXiv:hep-lat/0109019 



54 



[15] Ali Mostafazadeh, Supersymmetry and the Atiyah-Singer Index Theorem, J.Math.Phys. 



35 (1994) 1095-1124, [a£Xiv:hep-th/9309059l J.Math.Phys. 35 (1994) 1125-1138, 
arXiv:hep-th / 9309061] Supersymmetry, Path Integration, and the Atiyah-Singer Index 



Theorem, |arXiv:hep-th/9405048 [ 

[16 



[17: 

[18 
[19 
[20 

[21 

[22; 
[23 

[24 

[25; 

[26 

[27; 

[28 

[29 

[30 



G. V. Dunne, Aspects of Chern-Simons theory, Les Houches Lectures 1998, 



|arXiv:hep-th/9902115 



Joshua L. Davis, Per Kraus, Akhil Shah, Gravity Dual of a Quantum Hall Plateau Transition, 

JHEP 0811:020,2008, larXiv:0809. 18761 

Y. Zhang, et al., Phys. Rev. Lett. 96, 136806 (2006), 

Z. Jiang, Y. Zhang, H. L. Stormer, P. Kim , Phys. Rev. Lett. 99, 106802 (2007), 

Dipankar Chakrabarti, Simon Hands, Antonio Rago, Topological Aspects of Fermions on a 

Honeycomb Lattice, JHEP 0906:060,2009, larXiv:0904.T3T0] 

Michael Creutz, Four dimensional graphene and chiral fermions, JHEP 0804: 017, 2008, 
larXiv:0712.T20n 

A.Borici, Phys. Rev. D78 (2008) 074504, [arXiv: 0712.4 40T]. 

P.F Bedaque, M.I Buchoff, B.C Tiburzi, A.Walker-Loud, Phys. Rev. D78 (2008) 
017502.[ a7Xiv:0804.1145 ]. 

P.F.Bedaque, M.I.Buchoff, B.C. Tiburzi, A.Walker-Loud, Phys. Lett. B662 (2008) 449, 
[arXiv:0801. 3361]. 

L.B Drissi, E.H Saidi, M. Bousmina, 4D Graphene, Phys.Rev.D84:014504,2011, 
larXiv:110 6.5222l 

EI Hassan Saidi, On Flavor Symmetry in Lattice Quantum Ghromodynamics, J. Math. Phys. 
53, 022302 (20121. laFXTv: 1203.60041 

L.B Drissi, E.H Saidi, Dirac Zero Modes in Hyperdiamond Model, Phys. Rev. D84:014509, 2011, 
larXiv:1103. 13161 

L.B Drissi, H. Mhamdi, E.H Saidi, Anomalous Quantum Hall Effect of Graphene in Back- 
ground Fields, JHEP 026, 1110 (2011), [a?Xrv:1106. 55781, DOI: 10.1007, 

L.B Drissi, E.H Saidi, M. Bousmina, Electronic Properties and Hidden Symmetries of 
Graphene, Nucl.Phys.B829:523-533,2010, larXiv: 1008.4470 ^. 

Lalla Btissam Drissi, El Hassan Saidi, Mosto Bousmina, Graphene and Gousin Systems, in 
Graphene Simulation" Edited by: J.R. Gong, InTech Publishing, Rijeka, Croatia, (2011), 
larXiv:1108.1748| 



55 



K. S. Novoselov, Z. Jiang, Y. Zhang, S. V. Morozov, H.L. Stormer, U. Zeitler, J. C. Maan, G. 
S. Boebinger, P.Kim and A. K. Geim 2007 Science 315 1379, 

Yafis Barlas, Kun Yang, A. H. MacDonald, Quantum Hall Effects in Graphene-Based Two- 
Dimensional Electron Systems, Nanotechnology 23 052001 (2012). larXiv:1110.10691 
D. R. Hofstadter, Phys. Rev. B 14, 2239-2249 (1976), 

D. J. Thouless, M. Kohmoto, M. P. Nightingale, M. Nijs, Quantized Hall Conductance in a 
Two-Dimensional Periodic Potential. Phys. Rev. Lett. 49, 405-408 (1982), 
Mahito Kohmoto, Topological invariant and the quantization of the Hall conductance, Ann. 
Phys. (N.Y.) 160 (1985) 343, 

Giuseppe De Nittis, Giovanni Landi, Topological aspects of generalized Harper operators, To 
appear in: "The Eight International Conference on Progress in Theoretical Physics", Mentouri 
University, Constantine, Algeria, October 2011; Conference proceedings of the AIP, edited by 
N. Mebarki and J. Mimouni, (ajXiv:1202.0902[ 

D. H. Adams, Phys. Rev. Lett. 104, 141602 (2010) |arXiv:0912.2850] . 
D. H. Adams, Phys. Lett. B 699:394-397,2011, [ arX iv: 1008.2833]. 

Michael Creutz, Taro Kimura, Tatsuhiro Misumi, Index Theorem and Overlap Formalism with 
Naive and Minimally Doubled Fermions, JHEP 1012:041,2010, larXiv:1011.076n 
Michael Creutz, Taro Kimura, Tatsuhiro Misumi, Aoki Phases in the Lattice Gross-Neveu 
Model with Flavored Mass terms, Phys.Rev.D83:094506,2011, larXiv:1101.42391 
David H. Adams, Index and overlap construction for staggered fermions, Proceedings contri- 
bution for 28th International Symposium on Lattice Field Theory, Lattice2010, June 14-19, 
2010, Villasimius, Italy, Journal-ref: PoS (Lattice 2010) 073, i arXiv:110 3.619ll 
Michael Creutz, Confinement, chiral symmetry, and the lattice. Acta Physica Slovaca 61, No.l, 
1-127 (20111. [aFxrv:1103.3304t 

A.H.Castro-Neto et al. Rev. Mod. Phys.81, 109 (2009), 

Michael Creutz, Minimal doubling and point splitting, PoS Lattice 2010: 078, 2010, 
larXiv:1009.3 T54| 

L.B Drissi, E.H Saidi, M. Bousmina, J. Math. Phys. 52, 022306 (2011), 

L. Del Debbio, L. Giusti, M. L"uscher, R. Petronzio and N. Tantalo, Stability of lattice QCD 



simulations and the thermodynamic limit, JHEP 0602, Oil (2006), arXiv:hep-lat/0512021 



[47] By using periodicity property of the honeycomb supercell, one can show that (x) = -j^Q- 



56 



